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ABSTRACT 
The concept, instrumentation, and analytical techniques for the 
investigation of attenuation of wind generated water surface waves in 
a small scale real environment by separation of generation and propaga- 
tion by time in the dying wind situation were developed, and an empirical 


attenuation function applicable to a wide range of frequencies proposed. 
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TABLE OF SYMBOLS AND ABBREVIATIONS 


Definition Units 
discrete wave amplitude in time domain mm. 
complex Fourier amplitude mm. 
complex conjugate of a mm. 
equivalent bandwidth, M Af see 7 


arbitrary constant 

Computer Corporation 

diameter in. 

a discrete frequency sec 
aliasing or Nyquist folding frequency sec 
frequency sec 
feet 

acceleration of gravity ft-sec 
attenuation function 

Galvanized 


IBM computer library subroutine for fast 
Fourier transform 


an error function 

Hertz, cycles per second sec 

V -1 

inches 

International Business Machines 

frequency index 

angular spreading function 

attenuation constant sec? /nm 
time index 


pound, unit of force 


* (M) 


m, m(f), aE 


Mfd 


degrees magnetic 

attenuation coefficient 

mean water level 

microfarad, unit of capacitance 
number of elementary bands in B 
millimeters 

number of sample points 

nautical miles 

ohms, unit of resistance 

an integer 

a predicted energy density spectrum 
ratio of circumference of circle to diameter 
P(f£,x)/ S(£,x) 

root mean square 

Scientific Data Systems 


continuous energy density spectrum at a 
propagation position 


continuous energy density spectrum at a 
generation position 


average energy density in bandwidth B 
seconds 

signal generator 

a discrete time 

sample interval 

wave propagation time at group velocity 
length of wave record 

range of possible directions through which 


waves in a fetch can propagate to a down- 
wind point of interest 
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EE 


mm 


mm 


mm 


mm 


sec 


sec 


sec 


sec 


or 


-sec 


=Sec 


-Sec 


-SeCc 


limits of @ 
variation 
group velocity 


separation of two points along axis of 
wave propagation 


inches 


diode 


transformer 


capacitor 


resistor 


alternating current source 


ground potential 
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degrees 


ft-sec. 
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IT, INTRODUCTION 


The concept of wind-driven ocean waves as spectra of harmonic 
components which follow the laws of linear hydrodynamics was introduced 
in 1948 by Barber and Ursell (1). There has since been considerable 
progress toward the mathematical and physical description of ocean wave 
processes, with the ultimate objective being the prediction with reason- 
able certainty of future wave conditions at any geographical point on 
the basis of worldwide synoptic observations of both wave and wind fields. 
A summary of the major contributions to knowledge of ocean waves was 
compiled in 1967 by Cartwright (2). 

While theoretical models of wave generation, dispersion and geo- 
metrical spreading have been fairly well confirmed empirically and 
are used as the current basis of wave forecasting, the attenuation of 
waves during propagation is less well understood; and an attenuation 
function is consequently not applied in many wave —forecasting techni- 
ques. 

Some classical theories of wave attenuation, namely air resistance 
and turbulent friction, were critically examined by Darbyshire (4) and 
found wanting as singular explanations. The theory of non-linear scat- 
tering proposed by Phillips (6) appears promising and has been reason- 
ably supported by Snodgrass, et. al. (11). Conceivably, the answer 
may lie in a combination of these processes. 

As has been the case in studies of other environmental processes, 
it appears that the phenomenon under consideration must be completely 
described empirically before the question of theory can be resolved. 

So far, environmental studies of attenuation have been confined to low 


frequency ocean swell, where the range of measurable frequencies is 
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included in only one of several logarithmic decades of possible surface 
wave frequencies. The classical theories and the generalized model 
developed later in this paper call for the attenuation function to be 
exponential in propagation distance and frequency dependent. When the 
generalized model is applied to ocean wave statistics, the data are 
scattered within such a small band of the frequency scale that func- 
tions supporting almost all proposed theories can be fitted equally well 
to the data. 

Accordingly, the primary objective of this research was the collec- 
tion of attenuation data in a small scale real environment where the 
measurable frequencies are considerably displaced from those normally 
measured in ocean swell in order to give better definition to the at- 
tenuation function. This requirement for collection of data gave rise 
to the second objective, the development and evaluation of a system for 


recording and analyzing waves in a small scale real environment. 
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II. CONCEPTS 


A. COMPARISON OF ENVIRONMENTS 

Ocean wave attenuation studies involve the measurement of the energy 
of the relatively low frequency components of swell propagated over long 
distances from the area of generation. It is known that these low 
frequencies required high winds of appreciable duration for their develop- 
ment, hence the source of the swell waves can be located with fair ac- 
curacy in relation to the distance from the measuring instrument. A 
knowledge of wind or pressure conditions in the fetch can provide a 
reasonable estimate of the energy of waves in the fetch. 

High frequency waves can be generated by the prevailing low veloc- 
ity winds which are normally found over most of the ocean surface. For 
this and other reasons of distance and instrumentation, it would be ex- 
tremely difficult to make meaningful measurements of these frequencies 
with sufficient precision in the ocean to determine attenuation. For 
purposes of this paper, high frequencies shall be defined as those fre- 
quencies approaching, but not including, frequencies in which surface 
tension is the predominant restoring force rather than gravity. 

A lake or non-tidal estuary of moderate depth provides an ideal 
facility where experimental control is possible in the comparative 
measurement of high frequency deep water waves, although the scale of 
wind circulation in relation to the size of the lake or estuary provides 
no geographical separation between areas of generation and non-genera- 
tive propagation. On this scale, however, the mean wind may change 
rapidly in speed and/or direction. In a case where the wind speed 


diminishes rapidly, some frequencies generated at the higher wind speed 
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would receive no more energy and would propagate as swell. This con- 
cept, that fetch and swell can be separated by time just as there is 
the geographical separation in the ocean, is a fundamental assumption 


of this paper. 


B. THE GENERALIZED ATTENUATION FUNCTION 

Let S(f,x) represent the energy density spectrum computed from a 
wave record at a downwind (propagation) position, and let S(f,0) represent 
the energy density spectrum computed from a wave record at an upwind 
(generation) position, where f is the component frequency and x is the 
separation between two points of measurement along the axis of maximum 


propagation (the direction of the mean wind during generation). Then, 


S(f,x) = S(f,0) k(x,0,f) G(f,x) (1) 


where G(f,x) is the attenuation function, and k(x,9,f) denotes what 
portion of S(f£,0) would reach the downwind position if no attenuation 
were involved. This is a function of geometry and the beam pattern in 
the generating area. In this paper the model of Pierson, Neumann and 


James (9) will be used, where for all f 





9 
2 ie hes 
k(x,9,f) = i cos” 9 dO (2) 


The angles 8, and oh are measured at the edges of the fetch from the 
direction of the axis described above to the direction of the downwind 
measurement point. 


Let P(f,x) be a prediction of the spectrum at the downwind position 


which has not allowed for attenuation. Then, 


P(f,x) = S(f,0) k(x,0,f) H(f,x) (3) 
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where H(f,x) is an error function accounting for possible statistical 

errors in measurement at the upwind and downwind positions, and in- 

cludes differences between the angular spreading model and the actual 
angular spreading phenomenon. The model has proved to be fairly representa- 
tive of actual geometrical spreading, and H(f,x) should be at most a very 
slowly changing function of x. The spacial variation of H is thus ignored. 


Divide (3) by (1) 


R(E,x) = pee = HEE (4) 
Differentiate 
R'(£,x) =—SS— = -H(£,x) ee) (5) 
G*(£,x) 
Substitute (4) in (5) 
R'(£,x) = -R(f,x) —G-¢f2x) (6) 


G(f£,x) 


Now, assume that R(f,x) can be represented by: 


R(£,x) =C ger (7) 


where m = m(f), and C is a constant. This assumption appears to be 
valid from empirical data in ocean swell studies over a wide range of x. 


Differentiating (7) and substituting (7) in the result, 


R'(£,x) = Cme™™ = m R(f£,x) (8) 





then substitute (8) in (6), and 





G'(f£,x) = -m G(f£,x) (9) 


Ly 


which is a simple differential equation with solution 


G(f,x) =G(f£,0) e"" (10) 
Clearly, at the upwind measuring position, the function Go (£,0) = i 
and 

G(f,x) =e" (11) 


The attenuation coefficient m is a function of frequency and environ- 
mental conditions. Empirical evidence indicates that it increases with 
frequency, and from physical considerations it should reasonably in- 
crease from following wind to opposing wind conditions. Other environ- 


mental conditions may also be important. 
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IfIl,. INSTRUMENTATION 


Two resistance wave meters, adapted from a design by Pearlman (5), 
were constructed and used to collect data for this paper. A schematic 
diagram of the circuit used in both meters is shown in Figure 1. 

Pearlman's paper demonstrated that this type of wave sensor is es- 
sentially linear over a range of wave heights of a few inches in static 
and vertical oscillation tests and has errors of less than one percent 
when tested in orbital motion. Static tests of the modified design used 
in this research confirmed the linear characteristics of the instruments. 

The fixed resistors in the wheatstone bridge were standard 5% toler- 
ance resistors, but were matched by measurement on a resistance bridge 
to have values within 0.5 ohm of the resistance shown. The variable 
resistor was used to null the bridge at an initial immersion in the 
field prior to recording data. 

The signal generators used were Hewlett Packard Model 2000 audio 
oscillators. The recorder was a 14 channel Ampex CP-100 tape recorder. 
As the signal generators and the tape recorder both have one grounded 
terminal, the isolation transformer was required to prevent grounding 
one leg of the wheatstone bridge. 

Alternating current was used to minimize electrolytic effects on 
the electrodes. The full wave rectifier shown in Figure 1 demodulated 
this carrier frequency, passing a direct current signal to the re- 
corder that varied in response to the waves. The signal generators 
were operated at 3 volts (RMS) at about 2000 Hz; however, when operat- 
ing two wave meters on one recorder, it was necessary to have about 150 
Hz difference in the output of the signal generators to prevent a beat 


frequency from interfering with the wave records. 
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A dual beam oscilloscope was used to monitor the input and output 
of the tape recorder to insure proper recording of the wave signal. The 1] 
Mfd capacitor used across the rectifier output matched the input impedance 
of the tape recorder and oscilloscope, but the use of other equipments 
with different input impedances would require a change in capacitance to 
provide a smooth rectifier output without attenuating wave frequencies. 

The Ampex CP-100 tape recorder was operated with frequency modula- 
tion recording and reproducing amplifiers for the wave record channels 
at a tape drive speed of 1-7/8 inches per second. A third channel was 
used with amplitude modulation (direct) amplifiers for a voice time 
record. 

Details of the sensing probes are shown in Figure 2. The frames 
were designed to provide adequate rigidity, yet to minimize turbulence 
and reflections. The lower horizontal supports were designed to be at 
a depth greater than one-half wavelength of the anticipated waves. The 
fishing line shown was used as an insulator and to minimize the amount 
of electrode wire immersed. If the electrodes had extended to the lower 
horizontal support, the wheatstone bridge could not have been brought 
to near-balance by the variable resistor. In practice, it was found 
that the bridge could be nulled at an electrode immersion of no more 
than 2 centimeters. 

The probe frames were attached to 1 inch diameter galvanized pipes 
driven into the lake bottom by means of the 1% inch diameter sleeves 
shown in Figure 2. The meter circuits were mounted on top of the frame 
and connected to the oscillators and recorder ashore with Alpha #2311 


low capacitance shielded microphone cable laid on the bottom. 
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IV. DATA COLLECTION 


Roberts Lake, 1500 feet from Monterey Bay in Seaside, California, 
was used for the collection of data in support of this paper. The lake, 
shown in Figure 3, is uniformly 5 to 5% feet deep, with a soft mud bottom. 
Dense reed growth in the water around the circumference minimizes wave 
reflections. 

The Monterey Bay area was considered suitable for this small scale 
attenuation study because its climatological characteristics are favor- 
able to sudden changes in the wind. In addition to unpredictable changes 
in the wind field, diurnal heating and cooling of the Salinas Valley 20 
miles inland results in a sea breeze during the day which dies suddenly 
at dusk. This general characteristic is modified or not apparent when 
the central California coast is under the influence of transient high or 
low pressure regions. In three weeks of daily attempts to collect data, 
only one day produced the desired rapidly dying wind. 

Two sets of data, identified as data sets A and B, were collected 
on 21 November 1968. The wind speeds, as measured with a Casella pre- 
cision anemometer mounted eight feet above the lake surface, are shown 
in Figures 4 and 5. The Casella anemometer provides a mean wind speed 
by recording the revolutions of a three cup rotor on a digital counter. 
The wind data in data set A, Figure 4, were recorded for one-minute aver- 
ages, and in data set B, Figure 5, for half-minute averages. 

Details of wave sensor placement are shown in Figure 3. For data 
set A, the sensing probes were on an axis of 230°(M). Although the 
actual distance between probes was 400 feet, the direction of the wind 
(from 300°(M)) reduced the effective propagation distance to 138 feet. 


For data set B, the probes were 120 feet apart on an axis of 290°(M). 
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Prior to collecting the wave data, both probes were immersed 
simultaneously in relatively still water near the edge of the lake at 
two levels of immersion and one minute recordings of the signals were 
placed on magnetic tape for calibration of the wave meters. The average 
values of the one minute calibration signals were later used to scale 


the wave records to actual amplitudes in millimeters. 
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V. ANALYSIS OF DATA 


ae DIGITIZING 

The calibration records and selected portions of the wave records 
were digitized with an SDS 9300-COMCOR 5000 hybrid computer system. 

The Ampex CP-100 tape recorder was patched to the COMCOR 5000 analog 
computer through operational amplifiers included in the logic board of 
the computer and a low pass fiiter constructed from logic board components. 
The amplifiers were adjusted to amplify the maximum signal on the magnetic 
tape record to slightly less than the maximum input voltage of the analog 
computer. 

A schematic diagram and plot of the measured characteristics of the 
low pass filter are shown in Figure 6. To minimize aliasing in later 
spectral calculations, 128 Hz was selected as a filter cut-off frequency. 
Although frequencies of that order were far beyond the range of expected 
surface wave frequencies, noise could possibly have existed in that range 
from either inherent equipment noise or beats between the wave meter 
oscillators. 

To provide a minimum of two samples per cycle up to the cut-off 
frequency, the record was sampled at a rate of 256 samples per second. 


The digital computer sampling program used was developed by Post (10). 


Be CONVERSION OF DATA FOR IBM 360 COMPUTER SYSTEM COMPATIBILITY 

The core capacity of the SDS 9300 computer used for digitizing the 
data was far too small for spectral analysis of the wave records (see 
next section). Therefore, use was made of the much larger core capacity 
of the available IBM 360 computer system. The two computers do not share 


common symbology and storage parameters; however, a program developed by 
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Post (10) for conversion of data from one system to the other was used 
to transfer and convert the digitized data from the tape record made on 
the SDS 9300 computer to a tape record compatible with the IBM 360 


computer system. 


G. CALCULATION OF ENERGY DENSITY SPECTRA 
In a particular finite segment of a wave record of length T let 
A(t, ) represent the time series of N discrete wave amplitudes which are 
spaced at intervals of time At. Then the time of each sample amplitude 
is 
epee acne, k = 1,...,N CLD) 
The record is transformed from the time domain to the frequency 


domain by the Fourier transform 


N-1l 
A cel -2ITijk/N +_ 
ACE) = yy A(t, ) e sJ= 0,...,N-1 (13) 
k = 0 
where = j LANE are the component frequencies 


ACE) are the complex amplitudes of the f. 


J 
and os ae 

In the computer, this transform is accomplished by the IBM library 
subroutine HARM, which uses the fast Fourier transform algorithm developed 
by Cooley and Tukey (3), and requires that N be an integral power of 2. 

ACE.) is a two-sided amplitude spectrum where N/2 corresponds to the 
aliasing or Nyquist folding frequency, fi. which is the highest frequency 
analyzed. If the wave record had component frequencies above fo the ampli- 
tudes of these frequencies would be aliased into the spectrum ACE). A 


Nyquist frequency of 128 Hz was selected on the basis of the analog filter 
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characteristics as discussed in section A. The Nyquist frequency deter- 


mines the minimum sampling interval, Ac, by the relation 


At = 1/2£. (14) 


The energy density spectrum, AOE) is given by 


S(f,) = 2 Ace) Ax(£.) (15) 
| oe 


ANE 


Ae On Ni 


where Ax(e) is the complex conjugate of Ace) and the factor 2 accounts 
for the inclusion of only half of the two-sided symmetrical amplitude 
spectrum. The ea represent the average energy per unit frequency in 
an elementary frequency band of width A\f centered on f.. 

The NOE of one realization of an ergodic random process is a set 
of statistics, each of which is considered to have the Chi-square distri- 
bution with two degrees of freedom. The number of degrees of freedom 
determines the confidence limits for a statistic, and a meaningful compari- 
son Of statistics requires that the range of confidence limits be suffi- 
ciently small. The relation between the number of degrees of freedom 
and the confidence limits for the Chi-square distribution is given by 
Pierson and Marks (8). 

An accepted method of increasing the confidence in the statistical 
estimate of energy density is to average a number, M, of elementary 
frequency bands into a wider equivalent band of width Be When this is 
done, each bandwidth B has 2M degrees of freedom, As the two degrees 
of freedom in each elementary bandwidth reflect two independent samples 
in the time domain, the record length, tT required for the analysis is 
then specified by: 


T= Nee = 2f£, ut / B, (16) 
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Obviously, the increased confidence is bought at the expense of defini- 
tion. 

The parameters used for analysis of the wave records collected on 
21 November 1968 are summarized as follows: 


Es = 64 seconds 


N = 16384 
At = 1/256 seconds 
£ = 128 Hz 
fe 
Ba —eOnonHz 
e 
M = 32 


For 64 degrees of freedom, the 95% confidence limits are 1.5 and 
0.7, which means that for an infinite number of realizations of this 
particular wave process, 95% of the realizations of the spectrum would 
fall between 1.5 S(f,) and 0.7 Sie 

In the discussion of the results of the spectral analyses in the 
following section, confidence limits have been omitted for clarity, but 
their existence must be recognized. In fact, the computed S(£ ;) are the 


most likely values. 


SZ 


VI. RESULTS AND IMPLICATIONS OF SPECTRAL CALCULATIONS 


j 


In the text and figures of the discussion that follows, the time 
(hours:minutes:seconds) assigned to a spectrum is the mid-time of the 
G4-second wave record from which the spectrum was derived. For the wind 


data subset of data sets A and B, refer to Figures 4 and 5, respectively. 


A. COMPARISON OF SIMULTANEOUS UPWIND AND DOWNWIND SPECTRA 

Figure / shows the energy density spectra computed from simultan- 
eous realizations at the upwind and downwind wave probe positions at 
11:55:56 and 11:57:00. This comparison implies that spectra from down- 
wind probe data may be underestimated by a factor of approximately two. 
The most obvious cause for this under-estimation is an error in calibra- 
tion$ however, it is possible that the difference may be an effect of 


the geometry of the lake. 


Ine ATTENUATION INDICATED BY DATA SET A 

Figure 8 shows the energy density spectra at the upwind probe posi- 
tion at 11:57:00 when the average wind at the eight foot level had been 
four knots or more for at least 30 minutes, and at the downwind probe 
position 112, 176, and 240 seconds later. During the four-minute inter- 
val represented, the wind diminished to 1.3 knots. 

In the development of the generalized attenuation function earlier 
in this paper, the true energy density spectra at the generation and 


propagation positions were related by: 


S(£,x) = S(£,0) k(x,0,f) G(£,x) (1) 


and the attenuation function was specified by: 


mx 


Ged) — 8 Gi ky, 
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If it is assumed that the energy density spectra, S(f£.,0) and 
J 
S(f£.,x), derived from measurements of the upwind and downwind probe 
J 


positions, respectively, are equivalent to the continuous spectra, then 


G(f,x) = : : =e (17) 


it 


which, as the attenuation coefficient, m = m(f), is sought, can be 


written: AG Ie) k(x 4G) 
in y S65 ge 
BSE > a So (18) 


A valid comparison of energy densities for the determination of 
m(f .) can only be made when the identical wave components are present 
at both measurement points. Therefcre, SRB ou) measured at the upwind 
probe can only be compared with S(f 55x) measured at the downwind probe 
after the time, t, required for the energy of the component frequency 
“A to propagate over the distance x, the effective distance between probes. 


This is given for gravity waves by 
t= x 7 ve =4Gxn7f/e (19) 


where ve is the group velocity of the = and g is the acceleration of 
gravity. 

Table 1 summarizes the data and computations for Oa) for data set 
A. As the ae were computed at specific times, the frequencies whose 
propagation times corresponded to the time intervals between upwind and 
downwind spectra were computed by the inverse of equation (19) and 
rounded to the nearest % Hz to correspond to the center frequency of the 
band represented by a unless the calculated frequency fell midway be- 
tween center frequencies, then rounding was done in both directions and 


dual calculations made. 
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The second calculation of ay under each frequency uses 2 S(f.,x) in 
place of OED) to show the effect of the possible error in the down- 


wind spectra mentioned in section A. 


TABLE 1 


SUMMARY OF STATISTICS AND ATTENUATION COEFFICIENTS FOR DATA SET A 


t x Se eel EE ek, seer 
Upwind Downwind <aieem k(x,9,f) 1 
time time sec Ets sec mm /Hz mm /Hz ft 

bi ?57 700 11358 252 112 138 2.0 Lod see 648 .00473 
- .000348 

iee57200 11:59:56 176 138 $3.0 2.24 .95 asile) .0199 

0149 

Ss PRS, a99 0349 -0301 

sO251 

mees7 70086 12701:00 «6240 6«=6138)—S «4. 5 Sty) 95 .00265 + .0377 

-0327 


Now, as m = m(f), assume that m(f£) can be expressed by the equation 


m(f£) = K ¢P (20) 
where K is an attenuation constant and p is an integer. 
Table 2 shows the values of K for the values of m(f) given in 


Table 1 for p = 2,3,4. 
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TABLE 2 


ATTENUATION CONSTANTS FOR DATA SET A 


f m(f) m(£) pee pig p=4 
K K K 
seeo- fos _— seco an see fee 526 aw 
2.0 -00473 28.4 7.1 Syaays) Lay 
- .000348 -2.08 -.52 -.26 -.13 
3.0 00199 119. 13.4 4.41 1.47 
00149 89.3 ie, Sea 1.10 
See, 20301 180. 14.7 4.18 La2O 
0291 150. 1262 S29 1.00 
4.5 00377 226 1g 2.46 ~945 
TOD a7 196 9.6 7 a 472 
C. COMPARISON WITH OCEAN SWELL DATA 


Some energy density spectra of recorded Pacific Ocean swell from 
fetches at known distances from the swell observation point were made 
available by Wickham (12). Values of m and K for these data have been 
computed in the manner prescribed for these same values in Tables 1 and 
2 above, and are displayed in Table 3. The data for Table 3 may include 
situations in which the winds in the propagation areas were across or 
opposed to the direction of propagation, where the data for Tables l 


and 2 were for only the following wind situation. 
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TABLE 3 


ATTENUATION COEFFICIENTS AND CONSTANTS FOR OCEAN WAVES 





f m(f) p= 2 p = 3 p=4 

K K K 

-l -l 2 3 

sec nm sec /nm sec /nm sec /nm 
- LO5 £00275 Pe Bs: DisuheS Dihiad 
-095 ~00224 248 2.62 27.5 
095 .00261 ~289 3.04 3220) 
-085 - 00266 ~ 369 4.34 58.6 
075 -00378 -670 8.95 120. 
-055 -00115 - 380 6.93 L25i% 


D. RESULTS OF DATA SET A AND TABLE 3 

1. The possible error of a factor of 2 between spectra as dis- 
cussed in section A. does not have a significant effect on the deter- 
mination of attenuation for data set A. With the exception of the values 
at 7 = 2.0 Hz, the attenuation coefficients, m(f) are all of the same 
order of magnitude, and with the exception of the negative value for the 
second estimate of K for 2.0 Hz, the attenuation constants for each power 
of frequency are of the same order of magnitude. 

2. If multiplying St x) by 2 had a small effect, then variations 
within the confidence limits for the spectral estimates will also have 
small effect. 

3. It is quite evident that if an attenuation function can be 
described over a wide range of frequencies by a simple relationship such 


as proposed by equation (20), then it must be a function of ee 


E. DATA SET B 

The wind in data set B died more slowly than the dramatic change 
which allowed the extended calculations for data set A. For those short 
time periods where the wind declined rapidly, eS) were computed and 
are displayed in Figures 9, 10, and 11. As each figure represents only 
one attenuation observation, the corresponding values of she) and K are 
shown on the figures. 

Less confidence is placed in the results of data set B, as the 
results are more erratic and the spectra seem to be influenced near the 
frequencies of interest by a low frequency disturbance that is not evi- 
dent in data set A. This low frequency could have been due to a surge 


set up in the lake by the gusting of the wind during its decline. 


F. ADDITIONAL OBSERVATIONS 

The persistence of the energy density at 2.0 Hz is interesting, 
especially in view of the negative values obtained when 2 Se is 
substituted for S(f 5 5%) for the calculation of m(f ) and K. These 
negative values indicate that the energy in that particular frequency 
may be increasing with distance. A possible cause for this is that 2 Hz 
may be a low enough frequency to suffer negligible attenuation in such 
_a short distance. 

The second peak that occurs in both data sets at 6 Hz is also note- 
worthy. Possible explanations for this are that a frequency of 6 Hz is 
beginning to be influenced by surface tension effects, or the effect may 
be a result of non-linear interactions as have been reported by Pierson 


(7), 
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FIGURE 9. ENERGY DENSITY SPECTRA — DATA SET B 
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FIGURE 10. ENERGY DENSITY SPECTRA—DATA SET B 
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The decision to establish the Nyquist frequency at 128 Hz was not 
unwarranted, as examination of the spectra beyond the water surface wave 
frequencies show occasions of significant energy, particularly around 


60 Hz. 
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VII. CONCLUSIONS AND RECOMMENDATIONS 


The instrumentation, recording, filtering, digitizing and analyti- 
cal techniques employed in support of this investigation were extremely 
effective in providing well defined energy density spectra from waves in 
a small scale environment, 

The dying wind situation gives a definable separation of wind wave 
generation and swell propagation for high frequency waves when the wind 
dies rapidly and smoothly. 

Strong evidence is provided in support of a contention that an 
empirical attenuation function can be applied to a wide range of fre- 
quencies, and it is proposed that this function might be 


é 
G(f,x) = 5 " 


where K is approximately 3 ye ee for following winds and may increase 
to as much as 8 ee for opposing winds. Future extensive investi- 
gation may show the frequency exponent to not be an integer, but for the 
current state of definition of the attenuation function this integral 
representation is satisfactory. 

It is recommended that further studies be undertaken to provide 
better definition of such an empirical attenuation function, with emphasis 
placed on intermediate frequencies between the high frequencies investi- 
gated in this paper and the lower frequencies more frequently studied in 


the ocean. 
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